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We demonstrate the importance of electron correlation effects in the hyperfine structure constants 
of many low-lying states in 210 Fr and 212 Fr. This is achieved by calculating the magnetic dipole 
and electric quadrupole hyperfine structure constants using the Dirac-Fock approximation, second 
order many-body perturbation theory and the coupled-cluster method in the singles and doubles 
approximation in the relativistic framework. By combining our recommended theoretical results with 
the corresponding experimental values, improved nuclear magnetic dipole and electric quadrupole 
moments of the above isotopes are determined. In the present work, it is observed that there are 
large discrepancies between the hyperfine structure constants of the 7D 5 / 2 state obtained from the 
experimental and theoretical studies whereas good agreements are found for the other D 5 / 2 states. 

Our estimated hyperfine constants for the 8 P, 6 D, 105 and 115 states could be very useful as 
benchmarks for the measurement of these quantities. 

PACS numbers: 21.10.Ky, 31.15.-p, 31.30.Gs, 32.10.Fn 


I. INTRODUCTION 

Francium (Fr) is the heaviest alkali atom in the peri¬ 
odic table possessing 87 electrons. Therefore, the prop¬ 
erties of this system are expected to exhibit moderately 
strong correlation effects, and their determination calls 
for using powerful many-body methods. Like other alkali 
atoms, Fr atom is also being considered for many impor¬ 
tant experimental studies, prominent among them being 
the measurement of the electric dipole moment (EDM) 
due to parity and time reversal symmetries 0,0 and par¬ 
ity nonconservation (PNC) effects due to neutral weak in¬ 
teraction 0,0 and the nuclear anapole moment owing 
to its relatively heavy size. Like the EDM and PNC in¬ 
teractions, the magnetic dipole hyperfine interaction has 
a fairly strong dependence on Z , the atomic number of 
the system [5[ as it involves an electron interacting with 
the nucleus. Thus, theoretical investigations of hyperfine 
structures are necessary for EDM and PNC studies to 
test the accuracies of the wave functions in the nuclear re¬ 
gion 0 . On the other hand, comparison of theoretical 
results from various approximations with the experimen¬ 
tal values can provide a comprehensive understanding of 
the passage of electron correlations from lower to higher 
levels of many-body theory. This knowledge is essential 
to validate the theoretical results when the experimen¬ 
tal values are unavailable. Attempts have been made to 
investigate trends in the correlation effects in the calcu¬ 
lations of hyperfine structure constants of the S— states 
using lower order many-body methods 00 , but such 
trends have not been demonstrated explicitly for states 
having higher angular momentum in Fr. 

Here, we intend to show the variation in the trends of 
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correlation effects in the evaluation of hyperfine structure 
constants of as many as 17 low-lying states in Fr consid¬ 
ering relativistic second order many-body perturbation 
theory (MBPT(2) method) and the coupled-cluster (CC) 
method at various levels of approximation using the ref¬ 
erence states obtained by the Dirac-Fock (DF) method 
with F w_1 potential. In the present work, we have un¬ 
dertaken theoretical studies of of the hyperfine structure 
constants of 210 Fr and 212 Fr isotopes for two reasons. 
First, 210 Fr has been proposed as one of the most suit¬ 
able Fr isotopes for both the PNC and EDM studies 0 0 . 
To draw meaningful conclusions and to be consistent in 
the findings from comparisons between the theoretical 
and experimental hyperfine structure constants of differ¬ 
ent states, it is necessary to consider the results for as 
many as states as possible for both the isotopes. The 
second, but the essential reason behind considering the 
above two isotopes is that experimental results of hyper¬ 
fine structure constants for only a few selective low-lying 
states of 210 Fr |12|-|T5l and some other excited states of 
212 Fr are available llgl. 


II. THEORY 

The Hamiltonian describing non-central form of hyper¬ 
fine interaction between the electrons and nucleus in an 
atomic system is expressed in terms of tensor operator 
products as 0 

H hf = Y, Mi fe) -Og, (1) 

k 

where and Oare the spherical tensor operators 

with rank k (> 0) in the nuclear and electronic coor¬ 
dinates respectively. Since these interaction strengths 
become much weaker with higher values of k, we con¬ 
sider only up to k = 2 for the present case. Due to the 
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coupling between the electronic (J) and nuclear (I) angu¬ 
lar momentums, \"/IJ; FM f ) states are the proper bases 
with total angular momentum F = J + I and correspond¬ 
ing azimuthal quantum number M F and 7 representing 
rest of the unspecified quantum numbers. 

The energy splitting due to first order correction in the 
atomic state \JM) because of the hyperfine interaction 
is given by 

W^j = ( 1 IJ-,FM F \Y J ^ n k) .0[ k j\ 1 IJ;FM F ) 

k 

!, C} 

k ^ J 

X(I\\MW\\I)(J\\0$\\J) (2) 

which after expanding up to multipoles k = 2 , we get 

wgr = Wgj+wft, (3) 


where W F j and W F 2 are the contributions due to the 
magnetic dipole (Ml) with k = 1 and electric quadrupole 
(E2) with k = 2 interactions respectively. Traditionally 
it is expressed as 


and 


W: 


E 2 
F,J 


w¥,j = \ A hfK (4) 

f A 2 + jK — /(/ + l)J(J + 1) 
hf 21(21 -l)J(2J -l) 1 ’ 


with the hyperfine structure constants defined as 


A h f 


( J \\°hf\\J) 

UnQi — , 

yj J(J + 1 )( 2 J + 1) 


( 6 ) 


and 


B hf = 2 Q[ 


2 J(2 J — 1) 


(2J + l)(2J + 2)(2J + 3) 
x(J\\0$\\J), 


]l/2 


(7) 


where /in is the nuclear Bohr magneton, gj = (j-. nr 
and Q are the nuclear magnetic and quadrupole moments 
respectively, and K = F(F + 1) — /(/ + 1) — J(J + 1). It 
appears from the above expressions that the quantities 
Ahf/gi and Bhf IQ are independent of nuclear factors, 
but these values can slightly vary among isotopes due 
to different nuclear potentials seen by electrons in these 
systems. 


III. METHOD OF CALCULATIONS 

In the Block’s equation formalism, atomic wave func¬ 
tion (|d/ n )) of any of the considered states with closed- 
shell configuration [ 6 p 6 ] and a respective valence orbital 
n in Fr can be expressed as [Ts| 


where fi n is known as the wave operator that generates 
virtual excitations from the reference state |<f> n ) to form 
all possible excited configurations to give rise the exact 
state |'k„). The choice of |<I> ra ) can be crucial in obtaining 
accurate results, but we construct it using V N ~ X poten¬ 
tial by expressing |$ n ) = aj, |4>o) with DF wave func¬ 
tion |$ 0 ) of the closed-core [ 6 p 6 ] obtained by consider¬ 
ing the Dirac-Coulomb (DC) Hamiltonian (Hue) and is 
the working reference state for all the considered atomic 
states. 

To account for the electron correlation effects through 
H„ in the evaluation of |\l/„), we proceed in two suc¬ 
cessive steps. In the first step, we consider correlation 
effects within the closed-shell configuration [ 6 p 6 ], follow¬ 
ing which we take into account correlation of the valence 
electron coupling with the core electrons to generate the 
virtual excitations. For this purpose, we express f l n as 

fin = Ho + fi“, (9) 

where Ho (independent of n) takes care of correlation ef¬ 
fects from the closed-core while fi“ is responsible for in¬ 
corporating correlation corrections involving the electron 
from the valence orbital. In the perturbative expansion, 
we write [l 8 | 

2 

H n = (io) 

k =0 

with (lii 1 corresponds to wave operator with presence of 
k number of Coulomb interaction operator in it. The 
hyperfine structure constants in this method are deter¬ 
mined by calculating the reduced matrix elements of the 
O^l and operators (below denoted simply by O) 
using the following expression 

(fr n |o|fr») = ($ n inton„ |<f> n ) 

= (n) 

For the MBPT(2) method, i + j < 2 in the above sum¬ 
mations. 

In the Fock space CC approach, we define fi 0 = e T 
and fi“ = e T S n that yields the form 

|tfn> =e r {l + s n }|$ n ), (12) 

where T and S n are the excitation operators involving 
core and core-valence electrons respectively. Since our 
reference state is |$o), we take S v in normal order with 
respect to |$o) which, here, has been denoted by {S^}. 
It is worthwhile to mention that T is also in normal order 
with respect to |<J> 0 ) by construction. From a practical 
point of view we only consider the singles and doubles ex¬ 
citations through the CC operators, known as the CCSD 
method in the literature [19|, 123 , by defining 


|'F„) = H„|$„), 


( 8 ) T = T 1 +T 2 and S n = S ln + S 2n . 


(13) 









3 


Even though the CC operators have been approximated, 
we are able to account many important contributions 
from the triples and quadruples through the non-linear 
terms of Eq. m■ In order to understand role of these 
contributions in Fr for accurate evaluation of its hyper- 
fine structure constants, we also determine contributions 
only considering the linear terms of Eq. m by approx¬ 
imating 

I’l'n) = {1 + T\ + T 2 + Sin + <S , 2n}| < I > ra)j (14) 

which we refer to as LCCSD method. One more remark 
on this approximation is that it resembles the form of 
the singles and doubles configuration interaction (CISD) 
method. 

The hyperfine structure constants are determined in 
the CC method as 

(ttnlOIttn) = ($»|{1 + Sj}e T 'Oe T {l + S n }\<i> n ) 
($ n |{l + Si}e rt e r {l + S4|$ n ) ' 

(15) 

The above expression has two non-truncative series 
e T+ Oe T and e T+ e T in the numerator and denominator 
respectively. To account for the contributions that are 
significant from these series, we have used the Wick’s 
general theorem to divide these terms into effective one- 
body, two-body and three-body terms [3 . The effective 
one-body terms are dominant owing to the single parti¬ 
cle form of the above O operators, they are computed 
self-consistently and stored as intermediate parts before 
contracting with the corresponding S n operators. While 
the effective two-body and three-body terms are com¬ 
puted directly, but we also reuse the effective one-body 
terms to construct the effective two-body and three-body 
terms. To give an illustrate example of the adopted com¬ 
putational procedure in our calculations, we outline the 
important steps below. In the first step, we calculate 
the effective one-body terms of e Tt e T and store them as 
hole-hole (H-H), particle-particle (P-P), and particle-hole 
(P-H) blocks. The P-H block is first evaluated consid¬ 
ering only linear terms and then stored as an interme¬ 
diate effective P-H block which then contracted further 
with an extra T £ operator to get the H-P block. Conse¬ 
quently, the P-H block is obtained from the H-P block 
by multiplying appropriate phase factors and the proce¬ 
dure is repeated till the self-consistent results with tol¬ 
erance size 10 -8 achieved. Following this, the H-H and 
P-P blocks are constructed considering the direct con¬ 
tractions among the and Tj operators and along with 
the contractions of the H-P and H-P blocks with the T\ 
operators. As a result, the H-H and P-P blocks still con¬ 
tain terms from infinite series through the H-P and P-H 
blocks. In a similar fashion we compute e Tt Oe T but by 
slightly modifying the above strategy. Here we make use 
of the effective one-body terms of e Tt e T for the construc¬ 
tion of the effective one-body diagrams of e Tl Oe T and 
special attention has been paid to avoid the repetition 


TABLE I: List of different parameters of bases used in the 
present calculations. 


s 

V 

d 

/ 

9 

Ni 40 

39 

38 

37 

36 

a 0 2.0 x 10 -8 

2.5 x 10 -8 

2.5 x 10 -8 

2.1 x 10" 1 

2.1 x 10" 7 

/3 5.06 

5.04 

5.06 

5.08 

5.15 


of any diagram through the iterative procedure. We also 
replaced the Xj and its conjugate operators appearing in 
the effective two-body and three-body terms of e Tt Oe T 
by the effective P-H block and H-P block of e Tt e T to 
improve the results with the contributions from higher 
order CC terms. 

We also estimate contributions from the important 
triple excitations considering them perturbatively in the 
above property evaluation expression by defining a triple 
excitation operator 


qpert 
D 3 n 



ab,pqr 


(- H DC T 2 + H DC S 2n )Zn 


(16) 


where a, b and p, q , r indices represent for core and vir¬ 
tual orbitals respectively. Particles appearing in the sub¬ 
scripts are annihilated while those appear in the super¬ 
scripts are created in the course of defining excitation 
processes. We refer to this approach as the CCSD t 3 
method in this work. 

We use a recently developed basis function having the 
analytic exponential form with quadratic exponents to 
express the single particle wave functions. These func¬ 
tions for the orbitals having orbital angular momentum 
value l are given as 

Nl 

1 4>\r)) = r l J2 c l e ~ o “' ri \x{0,ip)), (17) 

y=i 


where |\(0, <p)) represents the angular momentum part, 
Ni corresponds to the total number of analytic functions 
considered in the calculations and a„ is an arbitrary co¬ 
efficient constructed satisfying the even tempering con¬ 
dition between two parameters ao and /? as 

u v = ao/3" -1 . (18) 


The motive for using the above analytic function is to de¬ 
scribe the atomic wave functions more accurately in the 
nuclear region of a heavy atom like Fr than the Slater and 
Gaussian type orbitals. We give the list of parameters in 
Table U used in the present calculations. 


IV. RESULTS AND DISCUSSION 

In order to demonstrate the propagation of electron 
correlation effects from the lower to higher orders of 
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TABLE II: Results for Ahf /gi and Bhf /Q (both are in MHz) obtained using the DF, MBPT(2), LCCSD, CCSD and CCSDt 3 
methods. Previously reported calculations using varieties of many-body methods are also quoted as “Others”. 





Ahf/gi 

values 



Bhf/Q values 

State 

DF 

MBPT(2) 

LCCSD 

CCSD 

CCSD t3 

Others 

DF 

MBPT(2) 

LCCSD 

CCSD 

CCSD t3 Others 

7 s 2 S 1 /2 

6531.06 

10186.68 

11457.29 

9916.20 

9885.24 

9927.69“ 
9947.07 6 
10155.40 c 
8432.65 d 






7 p 2 Pi/ 2 

707.70 

1169.41 

1522.50 

1281.84 

1279.56 

1308.67 6 












1265.77“ 












876.12“ 






7 P 2 P3/2 

55.68 

97.20 

124.80 

104.34 

104.28 

103.38 6 

118.03 

234.51 

300.47 

258.99 

259.73 231“ 







115.09“ 












96.28“ 






6 d 2 D3/2 

36.24 

44.10 

123.00 

85.02 

85.20 


26.19 

73.47 

100.17 

96.33 

98.61 

6 d 2 D 5/2 

14.28 

-23.04 

-84.36 

-58.92 

-58.82 


30.34 

97.09 

134.01 

129.38 

131.40 

8 s 2 Si /2 

1674.12 

2234.40 

2337.92 

2159.16 

2151.90 

2160.98“ 

2165.88 6 
2218.47“ 
2155.27 d 






8p 2 P\/2 

251.76 

381.06 

457.38 

402.84 

402.01 

408.68 6 












409.46“ 






8p 2 P3/2 

20.40 

33.96 

40.62 

35.04 

35.03 

34.25 6 

43.29 

78.74 

94.73 

83.81 

84.01 







39.64“ 






7 d 2 D3/2 

18.30 

16.14 

36.45 

30.84 

30.90 


13.23 

24.99 

29.28 

29.97 

30.53 

7 d “-D 5/2 

6.90 

-12.24 

-18.62 

-16.32 

-16.28 


14.57 

32.84 

38.56 

39.45 

39.92 

9s "S 1/2 

687.54 

881.28 

910.01 

852.78 

849.72 

852.39“ 

839.40 d 






10 s 2 S 1/2 

382.38 

468.12 

449.63 

455.82 

454.26 

402.97 d 






8 d "D 3/2 

9.38 

7.75 

16.62 

14.45 

14.47 


6.79 

9.47 

13.15 

13.43 

13.66 

8 d "-D 5/2 

3.46 

-5.34 

-7.62 

-6.90 

-6.88 


7.35 

12.62 

17.26 

17.59 

17.78 

9 d 2 D3/2 

5.08 

4.22 

9.13 

7.84 

7.86 


3.68 

5.01 

7.08 

7.12 

7.24 

9 d 2 D 5/2 

1.87 

-3.02 

-3.89 

-3.53 

-3.52 


3.96 

6.74 

9.21 

9.25 

9.35 

11 s 2 S\/2 

186.13 

238.72 

252.16 

231.13 

230.12 

208.ll d 







a [l5|, 6 [2l|, “[22|], d [ll|], e \2$. 


many-body methods, we present results using the DF, 
MBPT(2), LCCSD, CCSD and CCSD t3 methods for 
Ahf/gi and Bhf/Q of the first 17 low-lying states of 
210 Fr in Table [H] We find changes beyond third deci¬ 
mal places in these results for all the states except for 
the Ahf/gi values of the s-states, which are reduced by 
-1.141, -0.242, -0.095, -0.047 and -0.028 (in MHz) in 
the 7 S, 8 S, 9 S, 10S 1 and IIS' states, respectively, of 212 Fr 
compared to the results of 210 Fr. We account for these 
differences later while considering the results for 212 Fr. 
It can be seen from the above table that the MBPT(2) 
results are larger than the corresponding DF results for 


both the properties except for the Ahf/gi values of the 
D 5 /2 states. In fact, there are sign differences between 
the DF and MBPT(2) results of Ahf/gi in the D 5 / 2 
states implying that correlation effects are very strong 
for these states. It is found that the LCCSD method 
yields much larger values for both Ahf/gi and Bhf/Q 
compared to the other methods. The non-linear terms in 
the CCSD methods contribute substantially but with op¬ 
posite sign to reduce the LCCSD values. Triples through 
the CCSD t3 method also further reduce the Ahf/gi val¬ 
ues while they lead to small increases in the values of 
Bhf/Q . Thus, we conclude that the correlations effects 
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TABLE III: Contributions from individual LCCSD and CCSD terms to the Ahf/gi calculations in the 7S, 7 P, 6 D and 8 S 
states (in MHz). Extra and norm correspond to contributions from the non-linear terms of the CCSD method that are not 
mentioned explicitly and corrections due to the normalization of the wave functions respectively, cc means complex conjugate. 


CC term 

7 s 2 

LCCSD 

S1/2 

CCSD 

7 V 2 
LCCSD 

Pl /2 

CCSD 

7 P 2 
LCCSD 

P3/2 

CCSD 

6 d 2 D 3/2 
LCCSD CCSD 

6 d 2 

LCCSD 

D 5/2 

CCSD 

8 s 2 
LCCSD 

Sl / 2 

CCSD 

OTi + cc 

-440.401 

-437.580 

-11.058 

-18.540 

-1.050 

-1.740 

0.150 

0.042 

0.108 

0.066 

-101.161 

-101.52 

T/OTx 

7.422 

7.326 

0.042 

0.120 

0.006 

0.014 

- 0 

~ 0 

~ 0 

~ 0 

1.524 

1.536 

t/ot 2 

-17.701 

-19.321 

-0.486 

-1.260 

0.010 

0.018 

- 0 

- 0.002 

- 0.002 

0.002 

-4.092 

-4.230 

T/OT 2 

-72.300 

-75.721 

-4.548 

-5.628 

1.848 

1.920 

1.230 

1.068 

2.202 

2.226 

-16.356 

-16.548 

OSi v + cc 

3684.661 

2821.14 

486.901 

363.721 

31.956 

23.388 

57.438 

36.888 

18.901 

12.198 

494.341 

389.70 

OS 2v + cc 

1363.981 

1059.001 

225.720 

175.620 

28.698 

20.100 

-5.064 

-4.518 

-96.420 

-70.021 

308.040 

246.00 

T/OS 2v +cc 

-36.96 

-42.801 

-3.150 

-3.498 

-0.204 

-0.228 

-0.396 

-0.384 

-0.144 

-0.138 

-4.320 

-10.098 

T/OS 2v + cc 

-227.890 

-191.930 

-0.180 

-0.174 

0.099 

-0.060 

0.228 

0.126 

0.168 

0.036 

-47.240 

-40.91 

s\ v os lv 

519.661 

304.620 

83.821 

46.758 

4.590 

2.460 

22.962 

9.462 

6.300 

2.544 

36.486 

22.680 

S\ v OS 2v + cc 

288.301 

178.501 

53.310 

32.268 

5.412 

2.976 

3.498 

1.014 

-47.280 

-23.706 

18.186 

13.518 

s\ v os 2v 

353.880 

225.481 

31.326 

21.690 

0.570 

0.618 

23.358 

13.482 

8.724 

4.476 

74.580 

49.398 

Extra 


-90.476 


-12.198 


0.522 


- 2.022 


-2.208 


-56.382 

norm 

-495.831 

-267.311 

-46.836 

-24.738 

-2.604 

-1.326 

-16.596 

-6.378 

8.778 

3.072 

-91.911 

-46.723 


represented by the non-linear terms that correspond to 
triple and quadruple excitations are crucial for the accu¬ 
rate evaluation of the hyperfine structure constants in Fr 
atom. 

We also mention about previously reported calcula¬ 
tions of A h f /qi and B h f/Q in Fr in Table HD as “Oth¬ 
ers” [lH Il5 121H231 ]. Owusu et al. had investigated the 
core-polarization effects systematically in the S— states 
of 212 Fr using relativistic many-body perturbation the¬ 
ory equivalent to our MBPT(2) method. These val¬ 
ues are found to be smaller than our MBPT(2) results. 
SDpT method, analogous to our LCCSD approximation 
but with important triples effects in the wave function 
determination, has been employed in Refs. [HI, (2l| to 
evaluate the Ah / /gi values of a few low-lying states and 
are in close agreement with our CCSD and CCSD t3 re¬ 
sults. Dzuba et al had employed a restricted Hartree- 
Fock method in the relativistic framework and incorpo¬ 
rated correlation effects using many-body perturbation 
theory to investigate correlation effects in the hyperfine 
structure constants of a few low-lying states of 211 Fr 
[ 2 ^ . Heully and Martensson-Pendrill had employed a 
relativistic many-body perturbation method treating po¬ 
larization effects to all orders. Their results differ [23 
significantly from our calculations. 

Having demonstrated that substantial correlation ef¬ 
fects arise through the non-linear terms of the CCSD 
method above, it would be interesting to know whether 
these terms are more important for the accurate determi¬ 
nation of atomic wave functions which are necessary for 
PNC and EDM studies in Fr atom or the non-linear terms 
appearing in Eq. (fl5l) . In case of the second possibility, it 
would be advisable to avoid large scale computations in 


the evaluation of the amplitudes of the CC excitation op¬ 
erators. For this reason, we compare contributions aris¬ 
ing term by term of Eq. IT5l) in the LCCSD and CCSD 
methods. These values for Af,f / gi are given in Tables Hill 
m and 0 for each considered state of 210 Fr atom. Simi¬ 
larly, the calculated Bhf/Q values are given in Tables EO 
and IVIll Comparing contributions from individual terms 
of LCCSD and CCSD methods, we find the CC ampli¬ 
tudes are substantially changed by the correlation effects 
through the non-linear terms in the CCSD method. In 
fact, changes in the results due to normalization of the 
wave functions in both the methods (given as norm in the 
above tables) are also quite large. Contributions given 
as Extra are from the non-linear terms appearing in Eq. 
(fT5l) and found to be comparatively small. Clearly, the 
consideration of non-linear terms of the CCSD method in 
the calculations of atomic wave functions are very cru¬ 
cial, and they should be considered for obtaining high 
precision results and are very important for calculating 
quantities related to EDM and PNC studies in Fr. Even 
though Eq. © contains non-truncative series in the 
CCSD method, the linear terms contribute dominantly 
and are mainly responsible for determining the accura¬ 
cies of the final results. 

To understand the role of various correlation effects in 
the states belonging to different angular momenta, we 
discuss here the trends in the results through various CC 
terms for both the LCCSD and CCSD methods. It can 
be seen from these tables that maximum contribution to 
Ahf/gi come from OS± v and its complex conjugate (cc) 
terms, which correspond to all order pair-correlation ef¬ 
fects @,123, followed by OS 2 V + CC terms, which represent 
the all order core-polarization effects 6. [Hj . in the S, 
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TABLE IV: Contributions from individual LCCSD and CCSD terms to the Ahf/gi calculations in the 8 P, 7D, 9 S and lOS 1 
states (in MHz). Extra and norm correspond to contributions from the non-linear terms of the CCSD method that are not 
mentioned explicitly and corrections due to the normalization of the wave functions respectively. 


CC term 

8 P 2 
LCCSD 

Pl/2 

CCSD 

8 P 2 
LCCSD 

7*3/2 

CCSD 

7 d 2 
LCCSD 

7 I 3/2 

CCSD 

7 d 2 
LCCSD 

D 5/2 

CCSD 

9s 2 

LCCSD 

s 1/2 

CCSD 

10s 2 

LCCSD 

S 1/2 

CCSD 

OTi + cc 

-3.090 

-6.198 

-0.264 

-0.564 

0.084 

0.036 

0.054 

0.036 

-40.098 

-40.410 

-19.895 

-20.096 

Tl OT\ 

0.096 

0.036 

0.001 

0.001 

~ 0 

~ 0 

~ 0 

~ 0 

0.582 

0.594 

0.285 

0.291 

TloT 2 + cc 

0.138 

-0.276 

0.001 

0.003 

~ 0 

-0.001 

-0.001 

-0.001 

-1.626 

-1.680 

-0.807 

-0.833 

t}ot 2 

-1.764 

-0.252 

1.164 

0.942 

0.720 

0.540 

1.032 

0.972 

-7.476 

-7.524 

-3.846 

-3.872 

OSi v + cc 

121.861 

93.060 

8.670 

6.468 

11.700 

8.646 

3.606 

2.610 

162.780 

128.701 

74.151 

60.405 

OS 2v + cc 

71.940 

9.510 

9.666 

7.008 

1.524 

0.762 

-31.716 

-25.458 

122.341 

98.400 

60.806 

49.144 

Tl OS 2v + cc 

-1.074 

-1.212 

-0.072 

-0.084 

-0.150 

-0.156 

-0.126 

-0.054 

-3.504 

-4.080 

-1.761 

-2.051 

TlOS 2v +cc 

-0.540 

-0.060 

-0.036 

-0.018 

0.078 

0.048 

0.054 

0.030 

-18.321 

-15.962 

-9.020 

-7.874 

s\ v os lv 

14.760 

8.616 

0.924 

0.510 

1.920 

1.050 

0.486 

0.252 

9.630 

6.018 

3.964 

2.630 

S\ v OS 2v + cc 

8.880 

5.814 

0.948 

0.588 

2.160 

1.218 

-3.732 

-2.640 

3.342 

2.940 

0.941 

1.124 

s\ v os 2v 

9.540 

6.918 

0.306 

0.234 

6.084 

4.254 

2.526 

1.572 

29.508 

19.638 

14.666 

9.815 

Extra 


-6.408 


0.030 


3.366 


-1.584 


23.568 


-15.243 

norm 

-15.276 

-8.442 

-0.918 

-0.480 

-5.952 

-3.024 

2.250 

1.044 

-34.747 

-21.220 

-16.646 

-10.194 


TABLE V: Contributions from individual LCCSD and CCSD terms to the Ahf/gi calculations in the 115, 8 D, and 9 D states 
(in MHz). Extra and norm correspond to contributions from the non-linear terms of the CCSD method that are not mentioned 
explicitly and corrections due to the normalization of the wave functions respectively. 


CC term 

lls 

LCCSD 

2 s 1/2 

CCSD 

8 d 2 

LCCSD 

D 3/2 

CCSD 

8 d 

LCCSD 

2 d 5/2 

CCSD 

9 d 2 

LCCSD 

B 3/2 

CCSD 

9 d 2 

LCCSD 

D 5/2 

CCSD 

OTi + cc 

-10.583 

-10.704 

0.048 

0.023 

0.030 

0.021 

0.024 

0.013 

0.016 

0.011 

T/OTi 

0.150 

0.154 

- 0 

~ 0 

- 0 

~ 0 

~ 0 

- 0 

~ 0 

~ 0 

T}OT 2 + cc 

-0.429 

-0.443 

~ 0 

~ 0 

-0.001 

- 0 

- 0 

~ 0 

~ 0 

~ 0 

t ] 2 ot 2 

-2.111 

-2.113 

0.642 

0.235 

0.518 

0.476 

0.240 

0.180 

0.277 

0.251 

OS\ v + cc 

44.080 

35.231 

4.044 

3.009 

1.259 

0.907 

2.082 

1.510 

0.650 

0.458 

OS 2v + cc 

32.902 

26.513 

1.290 

0.831 

-14.290 

-11.668 

0.852 

0.570 

-7.452 

-6.095 

TlOSzv+cc 

-0.958 

-1.110 

-0.072 

-0.076 

-0.025 

-0.027 

-0.036 

-0.041 

-0.013 

-0.014 

T 2 OS 2v + cc 

-4.100 

-4.225 

0.006 

0.001 

0.014 

0.006 

0.004 

0.0002 

0.008 

0.003 

sl v os lv 

2.610 

1.667 

0.456 

0.253 

0.121 

0.063 

0.222 

0.118 

0.060 

0.030 

S\ v OS 2v + cc 

0.819 

0.776 

0.894 

0.553 

-0.635 

-0.527 

0.468 

0.292 

-0.273 

-0.220 

s\ v os 2v 

8.063 

5.353 

2.580 

1.855 

1.124 

0.725 

1.350 

0.967 

0.598 

0.386 

Extra 


-6.099 


-1.608 


-0.341 


-0.850 


-0.210 

norm 

-8.509 

-5.043 

-2.352 

-1.273 

0.795 

0.392 

-1.174 

-0.602 

0.366 

0.178 


P\/2 and D 3 / 2 states implying that the pair-correlation 
contributions are very important in all these cases. In the 
P 3 / 2 states, contributions from both these CC terms are 
almost same in magnitude while OS 2v + cc are the dom¬ 
inant contributors in the H 5/2 states. We also find that 
the trends of contributions in the evaluation of Bh f /Q are 
different from those of Ahf/gi- The core contributions 
coming through the excitation of the T operators are 
comparatively significant for the S— states. In general 


contributions through the single core excitations (gener¬ 
ated by T\) are found to be important. However, there 
are large core contributions through the T 2 operator, par¬ 
ticularly, for the ground state. 

We list the calculated and experimental values of Ahf 
and Bhf for 210 Fr and 212 Fr in Table IYTh! for which a 
few high precision measurements are available. Coc et al. 
had obtained the Ahf value for the ground 7 S state as 
7195.1(4) MHz using high-resolution spectroscopy (25|. 
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TABLE VI: Contributions from individual LCCSD and CCSD terms to the Bhf /Q calculations (in MHz) in the 7P 3 /2, 8 P 3 / 2 , 
6 D and ID states. Extra and norm correspond to contributions from the non-linear terms of the CCSD method that are not 
mentioned explicitly and corrections due to the normalization of the wave functions respectively. 


CC term 

7 p 2 

LCCSD 

P3/2 

CCSD 

6 d 2 

LCCSD 

D3/2 

CCSD 

6 d 2 

LCCSD 

D5/2 

CCSD 

8p 2 
LCCSD 

P3/2 

CCSD 

7 d 2 

LCCSD 

D3/2 

CCSD 

7 d 2 

LCCSD 

D5/2 

CCSD 

OTi + cc 

-2.237 

-3.690 

0.112 

0.033 

0.234 

0.141 

-0.564 

-1.202 

0.064 

0.028 

0.123 

0.081 

Tl OT\ 

0.010 

0.029 

- 0 

~ 0 

- 0 

~ 0 

0.002 

0.008 

- 0 

~ 0 

~ 0 

~ 0 

T/OTs 

-0.185 

-0.289 

-0.001 

-0.001 

- 0 

-0.001 

-0.051 

-0.094 

- 0 

~ 0 

~ 0 

~ 0 

T% OT2 

0.692 

0.888 

4.559 

4.332 

5.157 

4.926 

0.364 

0.248 

2.409 

2.060 

2.530 

2.162 

OSi v + cc 

67.760 

49.601 

41.611 

26.721 

40.121 

25.521 

18.390 

13.72 

8.487 

6.271 

7.665 

5.542 

OSxu + cc 

103.331 

88.490 

26.401 

39.060 

55.031 

68.220 

31.631 

27.34 

5.851 

10.56 

14.671 

19.77 

TlOSiv+CC 

-0.434 

-0.482 

-0.289 

-0.277 

-0.306 

-0.300 

-0.156 

-0.176 

-0.109 

-0.115 

-0.112 

-0.120 

Tl OS2v + cc 

-2.082 

-1.745 

-5.093 

-3.812 

-6.073 

-4.504 

-0.677 

-0.576 

-1.721 

-1.428 

-2.027 

-1.652 

s\ v os lv 

9.737 

5.216 

16.649 

6.858 

13.385 

5.409 

1.957 

1.090 

1.394 

0.760 

1.041 

0.543 

S{ v OS2v + cc 

14.081 

9.057 

0.240 

5.167 

8.190 

10.101 

2.534 

1.702 

2.830 

1.561 

3.293 

1.887 

S\ v OS2v 

-0.578 

-0.436 

3.294 

1.293 

1.873 

0.595 

0.153 

0.103 

1.620 

1.015 

1.472 

0.954 

Extra 


-1.786 


-0.006 


-0.045 


0.010 


0.010 


0.003 

norm 

-6.277 

-3.312 

-13.514 

-7.207 

-13.941 

-7.415 

-2.138 

-1.162 

-4.781 

-2.961 

-4.666 

-2.890 


TABLE VII: Contributions from individual LCCSD and CCSD terms to the Bhf/Q calculations (in MHz) in the 8 D and 9 D 
states. Extra and norm correspond to contributions from the non-linear terms of the CCSD method that are not mentioned 
explicitly and corrections due to the normalization of the wave functions respectively. 


CC term 

8 d 

LCCSD 

2 P > 3/2 

CCSD 

8 d 

LCCSD 

2 d 5/2 

CCSD 

M 

LCCSD 

2 P > 3/2 

CCSD 

9 d 

LCCSD 

2 d 5/ 2 

CCSD 

OTi + cc 

0.035 

0.017 

0.065 

0.044 

0.019 

0.010 

0.036 

0.024 

Tl OTi 

- 0 

- 0 

0.0001 

0.0001 

~ 0 

- 0 

0.0001 

~ 0 

tIot 2 

- 0.0001 

- 0.0002 

- 0.0001 

- 0.0002 

- 0.0001 

- 0.0001 

- 0 

- 0.0001 

T}OT2 

1.214 

0.992 

1.262 

1.035 

0.653 

0.523 

0.677 

0.545 

OS 3v + cc 

2.938 

2.185 

2.675 

1.926 

1.514 

1.096 

1.381 

0.972 

OS2v + cc 

2.429 

4.490 

6.338 

8.560 

1.226 

2.270 

3.246 

4.354 

TlOS2v+CC 

-0.051 

-0.055 

-0.053 

-0.057 

-0.028 

-0.029 

-0.028 

-0.030 

TlOS2v+CC 

-0.811 

-0.682 

-0.953 

-0.786 

-0.429 

-0.360 

-0.503 

-0.415 

s\ v os lv 

0.351 

0.183 

0.257 

0.133 

0.162 

0.085 

0.127 

0.063 

S\ v OS 2 v + cc 

1.350 

0.664 

1.364 

0.663 

0.756 

0.372 

0.748 

0.361 

slos 2v 

0.781 

0.517 

0.757 

0.515 

0.433 

0.287 

0.425 

0.291 

Extra 


-1.667 


-1.791 


-0.815 


-0.871 

norm 

-1.857 

-1.197 

-1.801 

-1.157 

-0.911 

-0.571 

-0.865 

-0.544 


In the same work, Coc et al had also given Ahf and 
Bhf values of the 7P 3 /2 state as 78.0(2) MHz and 51(4) 
MHz respectively. Using the two-photon excitation spec¬ 
troscopy of 210 Fr atom confined and cooled in a magneto¬ 
optical trap, Simsarian et al. had obtained the experi¬ 
mental Ahf value for the 8 S state as 1577.8(23) MHz at 
the ISOLDE facility [lj|. Fairly recently, Gomez et al. 
measured the Ahf value for the 9 S state as 622.25(36) 
MHz using a similar experimental technique [l5|. From 


another experiment, Coc. at al. have reported the Ahf 
value of the 7 Pi/ 2 state as 945.6(5.8) MHz [26|. A more 
precise value of Ahf in the 7P\/2 state has been obtained 
by Grossman et al. as 946.6(5.8) MHz [l2|. Soon af¬ 
ter that Grossman et al. extracted the Ahf values for 
the 7D 3 /2 and 7 D 5 / 2 states by measuring hyperfine split¬ 
tings as 22.3(5) MHz and —17.8(8) MHz respectively [14j. 
Here, they had neglected Bhf value of the 7 D 3 / 2 state 
while estimating a large Bhf value of 64(17) MHz for the 












7D 5/2 state. Similarly, for the 212 Fr isotope there are 
many experimental results available for the Ahf and Bhf 
values measured using various spectroscopic techniques 
and are presented in the same Table IVIIII For instance 
Coc et ah, adopting the same experimental technique as 
for 210 Fr, have measured the Ahf values of the 7 S and 
7 Pi/ 2i 3/2 states in 212 Fr [ 25 } 26]. It can be noticed from 
the tabulated results for the Ah / values that there is ex¬ 
cellent agreement between the present calculations with 
the experimental values of the 7P\/2,3/2 states whereas 
our results differ from the measurements by about 1.2% 
for the ground state. In another experiment Duong et al. 
had used stepwise laser excitation in collinear geometry 
with the on-line mass separator of the ISOLDE facility 
at CERN and measured the Aftrand Bhf values for the 
7 S, 7P and 8 P states of 212 Fr [27}. These experimental 
results are also in good agreements with our calculations 
for the above states. Arnold et al. had further extended 
this project to carry out measurements of the hyperfine 
structure constants of the IDS', IIS’, 8 D and 9 D states 
in 212 Fr [l6[> which are also given in Table rvTTTl It can 
be seen that these results and our calculations agree very 
well. 

As has been stated before, evaluation of the theoret¬ 
ical results for Ahf and Bhf require a knowledge of gi 
(i.e. gi and I) and Q values of the atom. Our calcu¬ 
lations using our CCSD t 3 method are the most rigor¬ 
ous theoretical results to date as they take into account 
more physical effects than previous calculations. The 
current best value of gi for 210 Fr was extracted by com¬ 
bining the experimental Ah / value of its 9 S state with the 
corresponding calculation using the SDpT method [l5} . 
However, if experimental results for hyperfine structure 
constants of any state are known to high precision, ex¬ 
traction of nuclear moments from these measurements 
can be justified. In reality, most of the measured Ahf 
values in 210 Fr are known quite precisely among which 
the ground state result is the most accurate (see Table 
IVIIII) . To infer gj value for 210 Fr from the Ahf results, 
we take the mean value from the data obtained com¬ 
bining the measurements with their corresponding calcu¬ 
lations using the CCSD t 3 method for all the states ex¬ 
cept for the D 5 /2 states. The reason for not considering 
results of these states is that the correlation effects in 
these cases are more than 100 % while for other states 
the principal contributions come from the DF values. In 
this approach, we obtain <?/ = 0.733765942 which cor¬ 
responds to fii = 4.40(5) of 210 Fr. This value is in ac¬ 
cordance with its earlier reported values as gj = 4.40(9) 
f28| and gi = 4.38(5) [k|. Unlike the case of Ahf, only 
two experimental values for Bhf in 210 Fr have been re¬ 
ported, among which the result for the 7 D 5 / 2 state might 
have been overestimated given that the wave functions 
of the D 5/2 states have an extremely small overlap in 
the nuclear region. Note that Bhf of the 7P 3 / 2 state 
is 51(4) MHz. Thus, combining the Bhf value of the 
7P 3 / 2 state with the corresponding calculation, we ob¬ 
tain Q = 0.196(15)6, where the uncertainty only from 


the measurement is taken into account, and the value 
estimated earlier was Q = 0.19(2)6 [28} for 210 Fr. The 
agreement between these two values is because the same 
experimental Bh / value has been used in both the results. 
By substituting these revised gi and Q values, we have 
evaluated the theoretical Ahf and Bhf values of 210 Fr 
and they are reported in Table IVIIII 

Keeping in mind the small differences among the cal¬ 
culations of the values of Ahf/gi and Bhf/Q between 
210 Fr and 212 Fr, as was mentioned earlier, we expect 
to observe the ratios between the Ahf and Bhf values 
from 210 Fr and 212 Fr for any given state to be almost 
equal to ratios of their gj and Q values, respectively, 
as per the demonstration in [29}. Considering all the 
experimental values known for the common states in 
both the isotopes, we find < 7 /( 212 Fr)/gj( 210 Fr) = 1.25(1) 
and Q( 212 Fr)/Q( 210 Fr) = —0.51(5). Experimental re¬ 
sults for Ahf and Bhf are reported for more states in 
212 Fr than 210 Fr. Excluding Ahf results for the D 5 / 2 
states owing to the reason stated previously, we ob¬ 
tain gi = 0.923070701 for 212 Fr when we combine its 
experimental Ahf values of the remaining states with 
their respective Ahf /gi calculations. This corresponds to 
/r/( 212 Fr) = 4.61(4). This, again, agrees with the previ¬ 
ously reported value gi( 212 Fr) = 4.62(9) [28[. In a similar 
procedure, we get Q( 212 Fr) = —0.10(1)6 which is same 
as that given in [28}. From these theoretical results, we 
get gj( 212 Fr)/g/( 210 Fr) = 1.26 and Q( 212 Fr)/Q( 210 Fr) a: 
—0.51, which are in reasonable agreement with the above 
estimated values from the measurements. 

Comparison of the theoretical and experimental Ahf 
and Bh / results quoted in Table IVIIII are quite satisfac¬ 
tory for almost all the states, but we find a large dis¬ 
crepancy between the experimental and theoretical Bhf 
values in the 7D 5 / 2 state. This result requires further 
theoretical and experimental verifications. Also, it was 
assumed that Bhf values of other D 3 / 2 5 / 2 states were 
negligible while extracting the experimental Ahf values 
of the corresponding states, however the present work 
shows that Bhf of the 7D 3 / 2 state is about 6 MHz. 


V. CONCLUSION 

In summary, we have employed many-body methods at 
different levels of approximation to study the magnetic 
dipole and electric quadrupole hyperfine structure con¬ 
stants of the first 17 states in 210 Fr and 212 Fr. This work 
demonstrates the importance of the inclusion of the non¬ 
linear terms in the coupled-cluster method, that accounts 
for the contributions from the triples and quadrupole ex¬ 
citations, for the accurate evaluation of the above quan¬ 
tities and for the studies of the violations of parity and 
time reversal symmetries in Fr. By combining the ex¬ 
perimental values with our corresponding calculations, 
we obtain gi = 4.40(5) and Q = 0.196(15)6 for 210 Fr 
and /ij = 4.64(4) and Q = 0.10(1)6 212 Fr. A reasonably 
good agreements between the theoretical and experimen- 
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TABLE VIII: Comparison between the theoretically determined and experimentally available Ahf and Bhf results (in MHz) 
of 210 Fr and 212 Fr. 




210 

Fr 



212 

Fr 


States 


Ahf 


Bhf 


Ahf 


Bhf 


Present 

Experiment 

Present 

Experiment 

Present 

Experiment 

Present 

Experiment 

7s 2 Si/ 2 

7253.45 

7195.1(4)“ 



9123.72 

9064.2(2)“ 
9064.4(1.5) 6 



7 V 2 Pi/2 

938.90 

945.6(5.8) c 



1181.12 

1189.1(4.6)“ 





946.3(2) d 




1187.1(6.8) 6 
1192.0(2) d 



7 P 2 P 3/2 

76.52 

78.0(2)“ 

50.91 

51(4)“ 

96.26 

97.2(1)“ 

97.2(1)*’ 

-25.97 

—26.0(2) 6 

6 d “D 3/2 

62.52 


19.33 


78.65 


-9.86 


6 d 2 D 5/2 

-43.16 


25.75 


-54.30 


-13.14 


8s 2 Si/ 2 

1578.99 

1577.8(23)“ 



1986.13 




8 P 2 P\/2 

294.97 




371.08 

373.0(1 ) 6 



8 p 2 P3/2 

25.68 


16.47 


32.34 

32.8(l) b 

-8.40 

—7.7(9 ) 6 

7 d 2 D3/2 

22.67 

22.3(5) / 

5.98 

Assume 0 ^ 

28.52 


-3.05 


7 d 2 D 5/2 

-11.95 

— 17.8(8)7 

7.82 

64(17)7 

-15.03 


-3.99 


9s 2 S 

623.50 

622.25(36) 9 



784.26 




10 s 2 Sl/2 

333.32 




419.27 

401(5)** 



8 d “D 3/2 

10.62 


2.68 


13.36 

13.0(6)'* 

-1.37 

Assume 0** 

8 d “-D 5/2 

-5.05 


3.49 


-6.35 

-7.1(6)'* 

-1.78 

- 2 ( 10 )** 

9 d “H 3/2 

5.77 


1.42 


7.26 

7.1(7)'* 

-0.72 

Assume 0** 

9d 2 D 5/2 

-2.59 


1.83 


-3.25 

-3.6(4)'* 

-0.94 

Assume 0** 

Us 2 S 1/2 

168.85 




212.39 

225(3)** 




“[2g], d [i3l, e 0, / B, ''Liiif- h M- 


tal results are obtained except for the electric quadrupole 
hyperfine structure constant of the 7D 5 / 2 state. Theo¬ 
retically predicted values for the hyperfine structure con¬ 
stants of many states including the 6 D states, in the 
present work, could be tested experimentally in the fu¬ 
ture. 


28,2013/4098. Computations were carried out using the 
3TFLOP HPC cluster at Physical Research Laboratory, 
Ahmedabad. 
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